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We investigate the interactions between two identical magnetic impurities substituted into a 
graphene super lattice. Using a first-principles approach, we calculate the electronic and magnetic 
properties for transition-metal substituted graphene systems with varying spatial separation. These 
calculations are compared for three different magnetic impurities, manganese, chromium, and vana¬ 
dium. We determine the electronic band structure, density of states, and Millikan populations 
(magnetic moment) for each atom, as well as calculate the exchange parameter between the two 
magnetic atoms as a function of spatial separation. We find that the presence of magnetic impurities 
establishes a distinct magnetic moment in the graphene lattice, where the interactions are highly 
dependent on the spatial and magnetic characteristic between the magnetic atoms and the carbon 
atoms, which leads to either ferromagnetic or antiferromagnetic behavior. Furthermore, through 
an analysis of the calculated exchange energies and partial density of states, it is determined that 
interactions between the magnetic atoms can be classified as an RKKY interaction. 
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Introduction 

For the last century, the vast increase in technological capabilities has been largely governed by the ability of devices 
to manipulate and control charge in various materials. However, charge is not the only property of electrons that 
can be exploited for technological use. Within the last decade, there has been a surge in the number of devices and 
materials that work to control, not only charge, but also the electron spin degrees of freedom. This evolution has 
brought about the area of spintronicsi"— , where devices made from spintronic materials are being investigated for 
applications ranging from quantum computing to enhanced nano-sized memory storage^”—, because they typically use 
less power and provide a dramatic increase in computer storage capacity and capability. 

The use for spintronic devices varies dramatically due to the multitude of interesting phenomena they exhibit. This 
includes spin relaxation and transport that may be used for quantum computation^”—. However, there are distinct 
challenges for spintronics, and a great need to completely understand the device applications for these materials^. 
Therefore, efforts in areas of experimental and theoretical physics, engineering, and chemistry are working to identify 
materials for spintronic application a^^d^ , as well as working toward understanding the fundamental interactions that 
govern these material propertiesi^. One large area of interest for these materials is that of two-dimensional (2D) 
materialsii^ii^, where fiat (or nearly-flat) organizations of atoms display extraordinary surface properties and may be 
useful for smaller and more accessible electronicsi^. 

The most popular example of a 2D material is that of graphen o^^J^ , which exhibits a large electronic mobility and 
thermal conduct ivit}^^”—!^ that provides the ability to conduct electrons. Recently, graphene has gained considerable 



FIG. 1: (Color Online) a) The 2D hexagonal structure of graphene and b) the calculated electronic band structure. 
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FIG. 2: (Color Online) a) The general substitution positions. 3-D view of the distorted graphene structure with two magnetic 
impurities separated by 1 carbon (b) and 6 carbons (c). 


attention in the scientific community due to the technological applications in electronic circuits and possible memory 
devicesi^ii^, which can allow for better transistors for faster processing and flexible circuits. 

Graphene consists of an individual layer of carbon atoms arranged in a standard honeycomb lattice (Fig. 

This configuration allows for graphene to have a large tensile strength and ability to conduct elec¬ 
trons through various p-orbital (tt) bonding24. Therefore, the tensile strength of graphene, relative to its weight 
and thickness, has been shown to be stronger than steel^^, which makes it a durable material for possible use in 
nanotechnological devices^. Furthermore, graphene has a specific crossover point in its electronic band structure, 
called the Dirac point, that occurs at the iF-point of the Brillouin zone (Fig. [Hb)), and provides the system a higher 
electron mobility. Because graphene is essentially a manifestation of surface states, it is considered to be a strong 2D 
topological Dirac material. 

Since graphene is being investigated as a base material for electronics, the current challenge is to establish it for use 
as a spintronic material. However, graphene is not magnetic by itsel f^^i^^ . Therefore, the ability to use graphene as a 
spintronic material is not very clear—. There have been a number of theoretical studies that have examined the viable 
of substituting magnetic atoms into graphene^”—. Recently, an experimental study showed that graphene deposited 
on an yttrium iron garnet thin film can induce a ferromagnetism in the graphene layer—. However, the mechanism 
that governs these interactions needs to be clarified. 

Another mechanism for possible magnetic induction in graphene is through direct magnetic substitution, where a 
carbon atom is removed from the graphene lattice and replaced with a magnetic atom. This can provide the necessary 
exchange energy and help clarify the underlying mechanisms. It is possible to provide direct substitution of magnetic 
atoms using scanning tunneling microscopy methods. The process of removing individual carbon atoms from the 
graphene lattice was recently demonstrated by Gomes et al4^. Furthermore, direct substitution of nitrogen and 
cobalt has been observe d^b^^ . Therefore, it is viable that magnetic atoms could be placed in the graphene lattice. 
This will be discussed further below. 

In this article, we present a comparative computational study examining the changes in the electronic and magnetic 
properties of transition-metal substituted graphene. Starting from a 128-atom supercell of graphene, two carbon atoms 
are replaced with either vanadium, chromium, or manganese impurities. The magnetic impurities are considered for 
various spatial separations of N carbon atoms (N = 1 - 6) in between the magnetic atoms along the zig-zag chain 
direction in graphene. Using density functional calculations with a full geometry optimization and spin polarization, 
we performed calculations to determine the total energy ground state, electronic band structure, density of states 
(DOS), and magnetic moment for each spin configuration. Overall, it is determined that the presence of magnetic 
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FIG. 3: (Color Online) a) Calculated exchange energy (change in total energy) for vanadium (black stars), chromium (blue 
squares), and manganese (red circles) substituted graphene, b) The partial density of states for the connector carbons (positive 
- spin up channel and negative - spin down channel), as well as that of pure graphene (dashed blue). Since this is for the five 
carbon case, there are three spin up and three spin down curves. This is because there are two pairs of carbons atoms that are 
spatially symmetry and therefore equivalent. 


atoms induces a magnetic state in the graphene, where the partial density of states shows a distinct increase in 
electron mobility at the Fermi level producing a metallic state. Additionally, we map the magnetic moment for each 
atom in superlattice to examine the magnetic exchange (or super-exchange) interaction characteristics throughout the 
graphene layer. Overall, this study provides clarity on how the magnetic atoms are communicating through graphene 
layer as well as helps in the understanding of the super-exchange interactions between the impurity atoms. 

Furthermore, through an analysis of the calculated exchange energy between the magnetic atoms, we find that 
the interaction alternates between a ferromagnetic and antiferromagnetic state as the magnetic atoms are separated 
spatially. This, in combination with the metallic state in the carbon, illustrates and suggests the presence of an RKKY 
exchange interaction. An RKKY (Ruderman-Kittel-Kasuya-Yosida) interaction is an indirect exchange coupling that 
is governed through interaction of localized spin moment and conduction electrons^”—. This produces two main 
criteria: 1) an alternating or oscillatory exchange strength dependent on distance, and 2) a distinct conduction 
of electrons typically in a metallic state. Therefore, this study provides avenues for the electronic and magnetic 
manipulation of graphene that can be utilized in the realization of spintronic and quantum computing devices. 


Results 

We conducted simulations on magnetically-substituted graphene lattice to see how the presence of magnetic atoms 
would affect the electronic and magnetic structure in the system. Therefore, we replaced two carbon atoms with 
vanadium, chromium, or manganese in a symmetric dimer configuration. Through DFT simulation of the formation 
energy, it was determined that graphene with magnetic atom substitution is energy favorable over a pure graphene 
state. To examine the formation energy, we determined the total energy for a pure graphene sheet with a unbound 
magnetic atom (total energy ^ -288 eV) and compared that to a graphene sheet with a bound magnetic atom and a 
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free carbon (total energy ^ -296 eV). Here, it is found that the magnetic substitution is energy favored by about 8 
eV. This is most likely due to the number of electrons that can be shared by the magnetic atoms over the carbon. 

To examine the electronic and magnetic interactions, we started with two atoms separated by a single carbon atom, 
and then separating them by one carbon along the zig-zag chain of graphene until there were six carbon atoms between 
the magnetic atoms (Fig. EJa)). This allows us to examine the spatial dependence of the magnetic interaction between 
the two impurities. In order to magnetically map the super cell, we calculate the Mulliken population over the entire 
lattice. Furthermore, we determine the electronic band structure and the density of states (DOS) for comparison to 
that of pure graphene (Fig. [TJa) and (b)). 


Impurity distortion 

The presence of the magnetic atoms produce a distinct distortion in the graphene lattice. This distortion is 
shown in Fig. EJb) and (c), where the three-dimensional view of the cell with 1 carbon and 6 carbon respectively. 
Through a geometry optimization, it is determined that the presence of magnetic atom impurities produces a dramatic 
distortion in the graphene lattice, which is consistent with previous first-principles calculations^^ and experimental 
measurements^!^. From this study, it can be determined that the bond lengths around the magnetic atoms are 
increased from the standard C-C length of 1.42 A to 1.87 A. This produces a distortion in the magnetic atom 
position and reduces the C-M-C bond angle from 120° to 92.9°. This is due to the destabilizing presence of the larger 
magnetic atoms. Further information on the distortion can be found in the supplementary material. 


Exchange energy and density of states 

Figure [3l^a) shows the exchange energy between the magnetic atom in each of the six configurations. The exchange 
energy is determined by examining the difference between the DFT calculated total energies for the FM and AFM 
states. According with the model, if the exchange energy is positive, then the interaction is FM. If the exchange 
energy is negative, then it is AFM. As shown in Fig. [Sj^a), all magnetic substitutions show a distinct RKKY variation. 
Although, only chromium and manganese fluctuate between FM to AFM, the vanadium interaction is constantly 
ferromagnetic and oscillates in the exchange strength. It should be noted that the magnetic ordering temperature is 
related to the strength of exchange coupling, where in this case, the ordering temperature of a spin dimer (FM or 
AFM) is dependent on the distance between the two magnetic atoms. 

Focusing mainly on chromium. Fig. shows the local density of states (LDOS) for the connector carbons in 

the chromium case (solid) as well as a pure graphene sheet (dashed). The DOS for all configurations is given in the 
supplementary material. Note that the density of states of pure graphene is zero at the Fermi level, while transition- 
metal substituted graphene has a distinct increase in the electron density of states at the Fermi level. This is a 
fundamental feature of Dirac materials where the presence of an impurity resonance at the Fermi level is produced as 
a gap forms in the Dirac cone^i.This provides a metallic DOS at the Fermi level that enables the conduction electrons 
to couple to the magnetic moment producing a super-exchange between the chromium atoms. This increase in the 
PD OS shows that the magnetic moments on the impurity atoms are influencing the electronic character of carbons 
between the magnetic atoms. Therefore, it is evident that the magnetic super-exchange is likely mediated through an 
RKKY coupling. 


Magnetic Induction 

Figure m^a) illustrates the absolute magnetic moment on each magnetic atom in their respective configurations. As 
shown, vanadium, chromium, and manganese has an average moment of around 1, 3.5, and 3 respectively. Here, 
/iB is the Bohr magneton. Figure SJ^b) shows the total magnetic moment for the entire supercell structure. Here, the 
vanadium case stays constant at about 2 /i^, which is mainly produced by the vanadium atoms maintaining a constant 
ferromagnetism. However, the net magnetic moment of the entire lattice structure tends to fluctuate for chromium 
and manganese as the magnetic atoms flip between FM and AFM ground states. Figure SJ^c) gives, for the example 
case of chromium, the magnetic moment on each of the connector carbons, where the inset details the pathway within 
the lattice. Therefore, there is an induced magnetic moment in the carbon, which means the simple super-exchange 
between the magnetic atoms is due to a complex network of interactions throughout the carbon atoms. 
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FIG. 4: (Color Online) a) Calculated magnetization for vanadium (black stars), chromium (blue squares), and manganese 
(red circles) substituted graphene, b) The calculated total magnetization for the magnetically substituted graphene, c) The 
magnetic moment for each atom between chromium atoms with the configuration for chromium substituted graphene separated 
by five carbon atoms. 


Electronic structure and magnetization mapping 

Figure [5] shows the electronic band structure for each ground state configuration of graphene with vanadium (a,d), 
chromium (b,e), and manganese (c,f) substituted atoms. The lower panels detail the lattice configuration and the 
magnetic moment for each atom in the super cell, where the magnetic atoms are denoted as the extreme cases. Overall, 
these panels illustrate the increased position of the magnetic atom in relation to each other by either 1,3, and 5 carbon 
atoms (Figure [5fa-c)) or 2, 4, and 6 atoms (Figure [5fd-f)). This figure illustrates the competition between the magnetic 
interactions between magnetic atoms and the carbon lattice. 

With respect to the electronic bandstructure, the competition between the interactions of the magnetic atoms 
and carbon atoms have distinct consequences on the presence of Dirac cone, which is the staple characteristic of 
pure graphene. It is clear that all of the odd carbon spacing (1, 3, and 5 carbon separations shown in Fig. [5fa-c)) 
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FIG. 5: (Color Online) The electronic band structure and configuration mapped magnetic moment in magnetically-substituted 
graphene with odd (1(a), 3(b), and 5(c)) and even (2(d), 4(e), and 6(f)) separations between two vanadium (a,d), chromium 
(b,e), and manganese (c,f) atoms. The magnetic moments of the magnetic atoms are not given by the scale. They have 
magnetic moments of 1, 3, and 3.5 jiB for vanadium, manganeses, and chromium, respectively. 


configurations seem to dissolve the Dirac cone. However, the Dirac cone seems to be restored within the spin down 
channels only when there exists an even carbon spacings (2, 4, and 6 carbon separations shown in Fig. Efd-f)) and 
a FM ground state between magnetic atoms. In the even cases for chromium (Figure [bfe)), the AFM character of 
the chromium-chromium interaction enables the whole sheet to be AFM. Therefore, the spin up bands are degenerate 
with the spin down bands. However, this highlights the competition between various interactions in the super lattice. 

The competition between M-C and C-C interactions seems to play a distinct role in the magnetic character for 
these configurations. In the vanadium cases, the magnetic atoms maintain an FM ground state throughout, which 
forces the C-C interactions in graphene lattice into either AFM and FM interactions. However, in the chromium 
cases, the C-C interactions between the magnetic atoms remain AFM, and the Cr-C interaction flips between FM 
and AFM. Therefore, the overall character of the Cr-Cr interaction changes and produces an RKKY-like fluctuation 
in the super-exchange. With manganese, the interactions are not trivial to determine. 
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Discussion 

As we examine magnetically-substituted graphene, we find that there are complex interactions between the elec¬ 
tronic, magnetic, and lattice degrees of freedom within the 2D structure. Within the electronic structure, the notable 
feature of a Dirac cone in pure graphene (shown in Fig. [I{b)) seems to disappear in most cases. However, when 
the magnetic atoms are separated by an even number of carbon atoms and interacting with FM character (as shown 
in Fig. Efd) and (f)), the Dirac cone seems to be restored within the spin down channel. From our calculations 
of the magnetic moment for each atom, it is clear that the carbon atoms within proximity to the magnetic atoms 
gain a distinct magnetization, as opposed to the zero magnetic moment of the pure graphene lattice. Therefore, 
the presence of a magnetic moment disrupts the electronic structure through a competition between valence (those 
used for magnetism) and conduction (those used for electronic transfer) electrons. However, when the carbon atoms, 
between the two magnetic atoms, have a FM character, that frustration seems to be released and the graphene lattice 
is able to restore the Dirac cone. We can determine that the Dirac cone is produced by carbon electron interactions 
because the Dirac cone is only restored in the spin down electron spectrum, which does not include the magnetic 
atoms. The Dirac cone details a distinct complexity in the competition between magnetic and electronic order in the 
superlattice. Although the significance of this phenomena is of technical nature, the present calculations do present a 
clear visualization of a Friedel-like spin density variation of induced moments on carbon atoms. We note that we have 
also calculated the case of non-magnetic titanium atom substitution, no spin polarization is observed in the system. 
Instead, a Friedel oscillation in the charge density is obtained. 

From Fig. [Ifb), we determine that the super-exchange between the two magnetic atoms is mitigated through both 
the spin up and spin down conduction electrons. This produces an increase in the PDOS at the Fermi level and 
explains why the carbon atoms are gaining magnetic moment when they are closer to the magnetic atoms (Fig. SJ^c)). 
Here, the sharing of electrons between the magnetic atoms and the lattice produces a net magnetic moment on the 
carbon sites. Based on the calculated Mulliken populations, we can infer that spin states for the magnetic atoms are 
Cr^+ {S = 3/2), Mn^+ (high-spin state S = 1), and {S = 1/2 or S' = 0) produced mainly from sharing three 

valance electrons with the graphene lattice. Vanadium is hard to gauge, since should have an S = 1 spin state. 
Therefore, we attribute its loss of magnetic moment to the possibility of a mixed valence state of and given 
a net magnetic moment of close to S = 1/2 (shown in Fig. SJ^a) and (b)). 

Therefore, the addition of magnetic atoms produces a net magnetic moment in the graphene lattice around the 
atoms, which leads to better conduction and an overall super-exchange. This exchange can be quantified by examining 
the change in the total energy of the FM (tt) state and the AFM (ti) state. Using a standard Heisenberg Hamiltonian 
for an independent dimer 


n = JSi-S2, ( 1 ) 

where Si is the quantum mechanical spin for sites i = 1,2, it can be determined that the change in energy {AE) from 
the FM and AFM state is simply proportional to the exchange energy J—^. Therefore, through an examination of 
the change in total energy of the two configurations, the super-exchange energy can be computationally estimated 
(shown in Fig. [3]^a)). This describes only the exchange between the two magnetic atoms. While there is an individual 
exchange between the induced magnetic moments on the carbon atoms, the interaction strength cannot be determined 
and can only be inferred as FM or AFM, based on the spin up or down nature of the magnetization. 

Depending on the coupling between the magnetic atoms and the carbons, magnetic frustration could be the cause 
of the variation between the FM and AFM characterization. As shown in Fig. [5fa)-(d), vanadium maintains FM 
character as the magnetic atoms are separated, which the local magnetic moments on the carbons in between orient 
themselves in response. This indicates that Jy-c (the exchange between vanadium and carbon) is dominant over Jc-c 
(the exchange between carbon atoms). Therefore, vanadium is able to force the carbon moments into a configuration 
that works to maintain Jy_ y >0. 

This effect is different for chromium (Fig. [Hb)-(e)) and manganese (Fig. [5fc)-(f)). In the chromium cases, Jc-c 
seem to dominate Jcr-Ci since the Jcr-Cr alternates through between FM and AFM, which is similar for manganese. 
These interactions seems to indicate the presence of an RKKY exchange through the graphene lattice. However, in 
order to have an RKKY interaction, the graphene needs to have a metallic DOS at the Fermi level, which is established 
in Fig. [3fb) by a distinct increase in the PDOS at the Fermi level for the carbon atoms between the magnetic atoms. 
Further work is being performed to try to characterize these interactions as well as the placement of other magnetic 
atoms. 

Overall, these calculations find a distinct RKKY exchange interaction. Since vanadium, chromium, and manganese 
are small transition-metal elements, the spin-orbit coupling is weak. The presence of spin-orbit coupling can produce 
a gap in the system forcing a standard ferromagnetic response. This was first investigated by Dugaev et al. using a 
standard exchange model with the spin-orbit coupling^, where they found that RKKY interactions should not be 
present due to this energy gap at the Fermi level. However, since we are using 3d transition-metal atoms, spin-orbit 
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coupling can be assumed to be negligible. Therefore, we find that the presence of magnetic transition-metal atoms 
produces multiple electron bands at the Fermi level beyond the standard Dirac cone crossing, and that there are 
distinct correlations between magnetic atoms (as shown in Fig. ED. 


Conclusion 

In this article, we present a comparative computational study examining the changes in the electronic and magnetic 
properties of magnetically-substituted graphene. Starting from a 128-atom supercell of graphene, two carbon atoms 
are then replaced with either vanadium, chromium, or manganese. The magnetic impurities are first considered with 
one carbon in between (1C) and separated until there are six carbons in between (6C) the magnetic atoms. Using a 
generalized-gradient approximation density functional theory with a full geometry optimization and spin polarization, 
we performed calculations to determine the total energy ground state, electronic band structure, density of states 
(DOS), and magnetic moment for each configuration. 

Through an analysis of the total energy for the different magnetic configurations, the super-exchange interaction 
between magnetic atoms is shown to alternate between antiferromagnetic and ferromagnetic state as a function 
of impurity distance. Additionally, we examine the spatially dependent magnetic moment for each atom in the 
superlattice and discover that the carbon atoms within close proximity to the magnetic atoms gain a distinct magnetic 
moment. The magnetic carbon atoms produce a complex network of exchange interactions. Through an analysis of 
the partial density of states for the carbon atoms, a distinct increase in electron mobility at the Fermi level is shown. 
Therefore, the calculation of a variable super-exchange with spatial distance and the transition to a more metallic state 
for the carbons in between the magnetic atoms indicates the presence of an RKKY-like super-exchange interaction 
through the graphene atoms. 

Further investigation is needed to understand the induced magnetic moment in carbon and the exchange network 
between the magnetic atoms. Given the advancements in scanning tunneling microscopy (STM) and nanoscopic en¬ 
hancements of atomic substitutions in graphene, experimental exploration of these interactions is a distinct possibility. 
The presence of variable magnetic exchange due to in-situ substitution and induced magnetism in graphene shows the 
increased possibility of using graphene for spintronic materials. These calculations are meant as the motivation for 
the use of magnetically-doped graphene and other 2D materials for the possible advancement for spintronic devices 
and quantum computation. 


Computational Methods 

Using Atomistix Toolkit^""—, we performed density functional theory (DFT) calculations within a spin-polarized 
generalized gradient approximation (SGGA) of PBE (Perdew, Burke, and Ernzerhof) functionalsi^ for a 2D hexagonal 
lattice of graphene supercell (Eig. [U^a)) with two magnetic impurities substituted such that they were separated by 
N carbon atoms through minimum distance, where N = 1-6 (Eig. [D^c)). The lattice consisted of a total of 128 
atoms and geometry was optimized using a force minimization method within an initial k-point sampling of 1x1x1. 
Eurther analysis was performed with a larger 3x3x1 k-point sampling once the geometry was minimized. The use of 
the smaller /^-sampling was for purposes of computational time given that extent and number of calculations being 
performed. Geometry optimizations performed with the higher /c-point sampling showed no significant difference. 

We implemented initial spin states of up-up (tt)? up-down (ti), and down-down (^^) configurations on the 
magnetically-substituted atoms and zero spin on the carbon atoms, which was established through non-substituted 
test runs on pure graphene. As expected, the configurations of tt and H were found to be degenerate. In all, we 
perform calculations on 54 total configurations (6 spatial configurations, 3 magnetic configurations, 3 different mag¬ 
netic atoms). The DET calculations were run using a tolerance of 0.3 /ieV with an energy cut off of 2 keV. As shown 
in Santos et al4^, the presence of a Hubbard U didn’t show any major differences for the smaller magnetic atoms 
of vanadium, chromium, and manganese. This is discussed in the supplementary materials, but is most likely due 
to the individual presence of the atoms. Therefore, since we are doing a comparative study, we did not implement 
an Hubbard onsite potential. It should be noted that while the presence of a Hubbard U would restrict some elec¬ 
tron movement by constraining the electronic orbitals, the general effects the trends for magnetization and electron 
conduction will still remain. 

In all, we calculated the total energy of each system to find the ground state configuration assuming collinear 
polarized spins and estimated the exchange parameters between the two magnetic atoms. Eurthermore, we determined 
the electronic band structure and density of states (DOS) to investigate the magnetic effect on the graphene Dirac 
cone at the K point in the Brillouin zone, and we also map the induced magnetic moments for each atom to examine 
the interactions throughout the graphene supercell. 
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